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ABSTRACT 

A semi-empirical, axisymmetric model of the solar minimum corona is devel- 
oped by solving the equations for conservation of mass and momentum with pre- 
scribed anisotropic temperature distributions. In the high-latitude regions, the 
proton temperature anisotropy is strong and the associated mirror force plays 
an important role in driving the fast solar wind; the critical point where the 
outflow velocity equals the parallel sound speed (v = cm) is reached already at 
1.5 R & from Sun center. The slow wind arises from a region with open field 
lines and weak anisotropy surrounding the equatorial streamer belt. The model 
parameters were chosen to reproduce the observed latitudinal extent of the equa- 
torial streamer in the corona and at large distance from the Sun. We find that 
the magnetic cusp of the closed-field streamer core lies at about 1.95 R Q . The 
transition from fast to slow wind is due to a decrease in temperature anisotropy 
combined with the non-monotonic behavior of the non-radial expansion factor in 
flow tubes that pass near the streamer cusp. In the slow wind, the plasma (3 is of 
order unity and the critical point lies at about 5 R Q , well beyond the magnetic 
cusp. The predicted outflow velocities are consistent with 5+ Doppler dimming 
measurements from UVCS / SOHO. We also find good agreement with polarized 
brightness (pB) measurements from LASCO/ SOHO and H I Lya images from 
UVCS/SOHO. 

Subject headings: Sun: corona — solar wind — Sun: UV radiation — MHD - 
Sun: magnetic fields 
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1. Introduction 

At the time of cycle minimum, the solar corona is more or less axisymmetric and stable 
for many months. The polar coronal holes are separated by an equatorial streamer belt that 
encircles the Sun. The fast solar wind originates from the coronal holes, while the slow wind 
originates from the vicinity of the streamer belt. The physical processes that drive the fast 
and slow wind are only partially understood. In particular, it is unclear why the Sun has 
such a bimodal outflow pattern with distinctly different physical conditions in the fast and 
slow wind. Does this distinction between fast and slow wind already exist at low heights in 
the corona, or does it arise only at larger distance from the Sun? What causes the transition 
from fast to slow wind as we approach the solar equator? What is the role of open and 
closed magnetic fields in this transition? To answer such questions, an empirical description 
of the corona is needed; i.e., a description of temperature, density, velocity and magnetic 
field as function of latitude, longitude and radial distance from the Sun. Such modeling is 
a necessary step in the study of the physical processes by which the fast and slow wind are 
generated. 

Axisymmetric, semi-empirical models of the solar corona have been developed by many 
authors (e.g., Pneuman & Kopp 1971; Yeh & Pneuman 1977; Steinolfson, Suess & Wu 1982; 
Washimi, Yoshino & Ogino 1987; Cuperman et al. 1990, 1993; Wang A. H. et al. 1993, 1998a; 
Lionello, Linker & Mikic 2001). For example, Sittler & Guhathakurta (1999) used empirically 
derived electron density profiles to construct a model of the magnetic field, outflow velocity, 
effective temperature, and effective heat flux. These parameters are derived by solving the 
equations for conservation of mass, momentum and energy, and the magnetic induction 
equation. The model provides an estimate of the large-scale surface magnetic field at the 
Sun, which is estimated to be 12-15 G. The authors predict that the large-scale surface field 
is dominated by an octupole term. More recently, three-dimensional models of the corona 
have also been developed (e.g., Linker, Van Hoven & Schnack 1990; Wang Y.-M. et al. 1997a; 
Riley, Linker & Mikic 2001). 

In this paper we develop a coronal model that synthesizes observational data obtained 
with instruments on the Solar and Heliospheric Observatory (SOHO) satellite, in particular 
data from the Ultraviolet Coronagraph Spectrometer (UVCS; Kohl et al. 1995). UVCS 
observations have shown that the minor ions in coronal holes have kinetic temperatures 
that are much larger than those for protons and electrons (Kohl et al. 1998; Cranmer et 
al. 1999). Here kinetic temperature refers to the total velocity dispersion of the particles, 
including thermal and non-thermal components. Detailed analysis of kinetic temperatures 
for different ions has shown that some ions have higher temperature than others, indicating 
preferential heating of certain ions. Furthermore, Doppler dimming analysis of spectral lines 
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such as O VI AA1032 and 1037 has shown that some minor ions have anisotropic velocity 
distributions: the velocity dispersion perpendicular to the (nearly radial) magnetic field in 
coronal holes is significantly larger than that parallel to the field. This high temperature 
and large temperature anisotropy of the ions is believed to be caused by dissipation of 
transverse waves (Dusenbery & Hollweg 1981; Isenberg & Hollweg 1983; Tu & Marsch 1997; 
Cranmer et al. 1999; Hu & Habbal 1999). A similar but smaller anisotropy may also exist 
for the protons (Cranmer et al. 1999). In the presence of a diverging magnetic field, charged 
particles experience an outward force that causes the perpendicular motion of the particles 
to be converted into parallel motion. Therefore, if the perpendicular energy of the protons is 
continually replenished by wave dissipation, the protons can maintain an anisotropic velocity 
distribution (T p ± > T p \\) and there is a net outward force on the protons. This so-called mirror 
force plays an important role in driving the fast solar wind (see review by Cranmer 2002). 

The purpose of the present paper is to include the effects of proton temperature anisotropy 
into a global model of the solar corona. The paper is organized as follows. In §2 we discuss 
observations of streamers and the relationship between the streamer and the slow solar wind. 
In §3 we present observational constraints on our global model, including temperature and 
density in the corona and magnetic flux at the coronal base. In §4 we present a method 
for solving the coronal force balance equations, taking into account the mirror force. In §5 
we present our model results for coronal magnetic structure and outflow velocity. In §6 we 
derive images of visible light polarization brightness and H I Lja intensities from our model, 
and we compare our results with observations. The main results of this work are discussed 
in §7. 



2. Streamer Structure and the Origin of the Slow Wind 

The K corona is produced by Thomson scattering of photospheric white light by free 
electrons in the corona. In a study of the solar minimum corona using the Large Angle and 
Spectrometric Coronagraph (LASCO) on SOHO, Wang Y.-M. et al. (1997a) found that the 
large-scale structure of the coronal streamer belt at 3 R & and beyond can be reproduced with 
a model in which the scattering electrons are concentrated around a single, warped current 
sheet that encircles the Sun. The angular width of this sheet is only a few degrees. The 
bright, narrow spikes seen in LASCO C2 and C3 images occur wherever the sheet is oriented 
edge-on in the plane of the sky (Wang Y.-M. et al. 1998b). In contrast, the latitudinal 
extent of the slow wind, as determined from in-situ measurements, is about ±20° (Suess et 
al. 1999), much larger than the thickness of the coronal plasma sheet observed with LASCO. 
This implies that the bulk of the slow wind originates from open field lines outside the 
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observed plasma sheet (Wang Y.-M. et al. 1998b). 

Sheeley et al. (1997) used time-lapse sequences of LASCO images to track the outward 
motion of small density enhancements ( "blobs" ) in the plasma sheet. These blobs originate in 
the high corona above the top of the helmet streamer, and slowly accelerate outward though 
the LASCO field of view (2.2-30 R Q ). Wang Y.-M. et al. (1998b) suggest that both the blobs 
and the plasma sheet represent closed- field material injected into the slow wind as a result 
of foot-point exchanges between the stretched helmet-streamer loops and neighboring open 
field lines. According to this model, the ejection of the blobs does not cause any permanent 
disruption of the helmet streamer, which remains in a stretched, quasi-equilibrium state with 
its cusp at 3 R & . 

Noci et al. (1997) and Raymond et al. (1997) observed the equatorial streamer with 
UVCS / SOHO at radii between 1.5 and 5 R & . Images of the streamer in O VI A1032 show 
two bright legs separated by a dark lane. This lane extends radially along the streamer axis 
up to 3 R Q (also see Strachan et al. 2002; Frazin, Cranmer & Kohl 2003). In contrast, H I 
Lya and white-light images do not show such a dark feature. This led Noci et al. (1997) to 
attribute the dark lane to a reduced oxygen abundance along the streamer axis. Raymond 
et al. (1997) further suggested that this abundance anomaly is due to gravitational settling 
of oxygen in the static proton/electron plasma of the closed-field streamer core. If this 
interpretation is correct, the closed field lines must extend up to 3 R & , consistent with the 
model of Wang Y.-M. et al. (1998b). 

The O VI images for July 1996 show that the streamer legs converge towards the equator 
in the range 2-2.5 R Q from Sun center. This suggests that the streamer cusp may be located 
at about 2 R & , somewhat less than implied by the model of Wang Y.-M. et al. (1998b). 
We propose the following scenario. The blobs observed by Sheeley et al. (1997) originate 
at the cusp (2 i?©), but become clearly visible only at somewhat larger height (3-4 R & ). 
The plasma within the blobs originates in the closed-field streamer core, which is affected 
by gravitational settling. Therefore, the blobs have low oxygen abundance compared to the 
surrounding slow solar wind. As the blobs move out into the region between 2 and 3 R & , 
they contribute to reduced O VI emission along the axis of the streamer. In this way, the 
effects of gravitational settling are transmitted to larger height via the blobs. In contrast, 
the open field lines farther away from the streamer axis always have a slow outflow and 
therefore are unaffected by settling. 

Strachan et al. (2002) measured 5+ outflow velocities in streamers using the Doppler 
dimming effect (also see Habbal et al. 1997). In a latitudinal scan at 2.33 R Q , Strachan et 
al. found no measurable outflow velocity within the streamer (v < 20 km/s), but a steep rise 
in outflow velocity occurs just beyond the bright streamer legs. The latitudinal width of the 
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streamer, as defined by the point where the outflow velocity equals 100 km/s, is about ±20° 
(see Fig. 4d of Strachan et al. 2002), similar to the observed width of the slow wind at large 
distance from the Sun (Suess et al. 1999). 

3. Observational Constraints on Global Model 

This section describes the observational constraints on coronal density, temperature, 
and magnetic flux to be used in §4. 

3.1. Electron Density 

The polarized brightness (pB) of the K corona can be used to measure the coronal 
electron density (van de Hulst 1950). Here we use the results of Guhathakurta, Holzer & 
MacQueen (1996) and Sittler & Guhathakurta (1999), who used Skylab data from 1973-1974 
to derive electron density profiles for coronal holes and streamers at cycle minimum (see 
also Newkirk 1967; Allen 1973; Munro & Jackson 1977; Saito et al. 1977). Their results 
are expressed in terms of radial profiles, one for the polar coronal hole, another for the 
equatorial streamer. They used the following expression for the electron density within each 
component: 

N e (r) = ai e a2/r r- 2 [l + a 3 /r + a A /r 2 + a 5 /r% (1) 

where r is the heliocentric distance in units of R Q , and the parameters a±, ■ ■ • , 05 are given by 
Sittler & Guhathakurta (1999). The density profiles are shown in Figure 1. The Skylab data 
span the height range 1 — 5-R Q and were further extrapolated using data from the Ulysses 
mission (Phillips et al. 1995). These measurements correspond to select days when the 
equatorial streamer belt was seen approximately end-on. This explains why their streamer 
densities are somewhat larger than those of Saito et al. (1977), who presented streamer 
densities averaged over many days. 

3.2. Electron Temperature 

One method for inferring the electron temperature T e is to use spectral line intensity 
ratios of two lines from the same ion. Wilhelm et al. (1998) used data from SUMER/ SO HO 
to estimate T e at radii 1.03-1.6 R Q in coronal holes at the last cycle minimum (1996-1997). 
They use Mg IX 706 A and 750 A, and they conclude that in both the plume and inter-plume 
regions the electrons barely reach the canonical temperature of 1 MK. Moreover, T e (r) falls 
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off rapidly with height. Although their temperature estimates depend on atomic data that, 
according to the authors, could be improved, their work suggests that in polar hole regions 
the electrons are significantly cooler than the ions (see Landi et al. 2001, for a quantitative 
analysis on uncertainties of T e estimates, derived from Be-like line ratios and using different 
theoretical methods). 

The rates of ionization and recombination of coronal ions decrease rapidly with distance 
from the Sun. Therefore, the charge state of the solar wind is determined in large part 
by the electron temperature T e in the inner corona where the ionization and recombination 
times are still short compared to the solar wind expansion time. In situ measurements of 
the solar wind charge state can be used to estimate the coronal electron temperature. Ko 
et al. (1997) derived polar T e profiles from observations with SWICS/Ulysses. Their results 
can be approximated as a combination of two power laws (see Cranmer et al. 1999): 

T e (r) = 10 6 (0.35 r 1A + 1.9 r~ 6 - 6 ) _1 [K]. (2) 

Raymond et al. (1997) and Li et al. (1998) studied the ionization balance of various ions in 
a streamer observed with UVCS/ SOHO in July 1996. The results indicate that T e reaches 
a maximum value of about 1.6 MK in the streamer core. Due to the high density and low 
outflow velocity in streamers (e.g., Frazin 2002), we expect that protons and electrons are in 
thermal equilibrium with each other, so the electrons may be used as a proxy for the protons. 
Raymond et al. (1997) and Li et al. (1998) show that their observations are compatible with 
hydrostatic equilibrium in the streamer core (Gibson et al. 1999, also see). In this paper we 
derive the electron temperature T e (r) from the observed electron density N e (r), assuming 
hydrostatic equilibrium. 

In this paper we use a generalization of expression (2) to describe the radial variation 
of electron temperature at the pole and the equator: 

Te(r) = T ° a + br° a +\l - b)r-f (3) 

where T is the temperature at the coronal base (r = 1). The values of the parameters a, b, a 
and (3 are given in Table 1, and the profiles are shown in Figure 2. The electron temperature 
model at the poles (dashed line) has been adjusted to approximately fit the observational 
data of Ko et al. (1997) and Cranmer et al. (1999) [triangles). The electrons are assumed 
to have a Maxwellian velocity distribution. The electron temperature at the equator (thick 
line) is assumed to be equal to the proton temperature, discussed in the next section. 
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3.3. Proton Temperature 

Ion temperatures can be estimated by measuring the width of coronal emission line 
profiles. Due to rapid charge exchange between protons and neutral hydrogen, the latter 
can be used as a proxy for the protons. Allen, Habbal & Hu (1998) studied the coupling 
between neutral hydrogen and protons by treating the hydrogen atoms as test particles in a 
proton-electron background. Their work indicates that the H I velocity distribution reflects 
the proton distribution at radii up to about 3 R & in the polar regions and even higher in the 
equatorial regions. Above 3 R Q in the polar regions, the H I velocity distribution "follows" 
the proton distribution and reaches temperatures about 20% higher than the protons (see 
also Olsen, Leer & Holzer 1994; Olsen & Leer 1996). These results strongly support the 
use of H I Lja profiles to measure the velocity spread of the protons along the line of sight 
(LOS). The observed velocity spread includes both thermal and non-thermal components 
(such as wave motions), and may also include a contribution from solar wind expansion in 
the direction along the LOS (most emission originates near the point of closest approach 
to the Sun, but there is some contribution from regions behind and in front of the plane 
of the sky where the solar wind velocity has a component along the LOS). Therefore, the 
velocity spread derived from the observed line width provides only an upper limit on the 
proton temperature: T p < (mn/kB)Vy c /2, where Vi/e is the observed 1/e velocity spread, 
van is the hydrogen mass and ks is Boltzmann's constant. At the time of cycle minimum, 
the polar coronal holes are very large, and polar observations at r > 1.5-R© are unaffected 
by low-latitude streamers. 

UVCS / SOHO observations of heavier ions such as 5+ show that the velocity distribu- 
tions of these ions in coronal holes are highly anisotropic: the velocity spread in the radial 
direction (as derived from Doppler dimming analysis) is much smaller than that in the tan- 
gential direction (Kohl et al. 1998; Cranmer et al. 1999). It is unclear whether this anisotropy 
also exists for the protons; Doppler dimming analysis of H I Lyman lines does not provide 
strong constraints on the proton parallel velocity in coronal holes. The 5+ anisotropy is 
believed to be due to the damping of transverse waves (e.g., due to ion-cyclotron resonance), 
and this perpendicular heating may also occur for the protons. In this paper we assume that 
the protons indeed have an anisotropic velocity distribution in the low-density polar regions 
(T p ± > T p \\), but not in the equatorial region (T p ± = T p »). Since the magnetic field over the 
pole is approximately radial and perpendicular to the LOS, the observed H I Lja line width 
provides an estimate for the proton perpendicular temperature. Furthermore, we assume 
that the proton parallel temperature equals the (isotropic) electron temperature, T p \\ = T e , 
so that the radial variation of T p \\ at the pole and equator is given by equation (3) with 
parameter values given in Table 1. The anisotropy in proton temperature over the pole is 
consistent with in-situ measurements, which show that T p > T e in the solar wind (Marsh et 
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al. 1982; Pilipp et al. 1987). 

Kohl et al. (1998) and Cranmer et al. (1999) analyzed H I Lya line profiles observed 
in polar coronal holes during the past cycle minimum (1996-1997). Their results indicate 
proton perpendicular temperatures up to 6 MK. The observed velocity width (Vi/ e ) increases 
rapidly with height from about 190 km/s (T p± ~ 2.2 MK) at 1.5 R Q , to 240 km/s (T p± ~ 3.5 
MK) at 2.5 R & , and then slowly decreases to 250 km/s {T p± ~ 3.8 MK) at 4 R Q . The results 
are consistent with earlier measurements from UVCS / Spartan by Kohl, Strachan & Gardner 
(1996), who found peak temperatures of 5-6 MK at 2.25 R Q and 3.5 MK at 3.5 R & in a polar 
coronal hole in 1993. From these several observations, measured temperatures at different 
heights are shown by asterisks in Figure 2. 

UVCS / SOHO observations in the equatorial regions indicate a roughly constant proton 
temperature within the core of the streamer belt. We reanalyzed UVCS data from a super- 
synoptic campaign during July 1996 (Raymond et al. 1997), and find Lya line widths of 
order 185 km/s (T p± ~ 2.0 MK) at 1.5 R Q and 195 km/s (T p± ~ 2.2 MK) at 2.6 R Q . At 
larger heights, the line widths decrease to 170 km/s (T p± ~ 1.75 MK) at 3 R Q and 150 
km/s (T p j_ ~ 1.35 MK) at 4.5 -R . The measurements are shown by the crosses in Figure 
2. Streamer observations by Kohl et al. (1997), obtained about a month after the Raymond 
observations, exhibit very similar values (2.2 MK at 2 R Q , 1.5 MK at 4 R Q ). 

In this paper we use the following expression for the proton perpendicular temperature: 

a + 1 (r - l) 2 e-i*- 1 )/^ 

Tp±(r) = T °a + br« + (l-b)r-P + Tl (r max - If e -(*w,-i)/Ar ' (4) 

The first term is similar to that used for the electron temperature, but the second term 
allows us to impose a further increase in temperature over a limited range of heights (of 
order Ar), consistent with UVCS observations. The values of the parameters at the pole 
and at the equator are given in Table 1. 

In summary, our temperature models present the following main features: (1) at the 
pole, strong anisotropy for protons, and electron temperature much lower than T p ±; (2) at 
the equator, hydrostatic equilibrium isotropic velocity distributions and thermal equilibrium 
between species. Note that T and (5 are the same for all temperature models, so that 
T p ± = T p \\ = T e at low heights in the corona. 



3.4. Magnetic Flux at Coronal Base 



Figure 3 shows the longitude-averaged radial magnetic field B r as function of latitude 
at the solar surface (r = 1). These data are derived from the Kitt Peak synoptic map for 
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July 1996, near the time of cycle minimum. There are no large active regions on the Sun at 
this time. The full curve in Fig. 3 is a fit of the form 

B r (6) =B cos p e, (5) 

where p = 7 and B = +10 G. Note that B r decreases monotonically from about +10 G 
in the North to -10 G in the South. At latitudes between ±50 degrees, the average radial 
field is very small (less than 2 G). Therefore, the large-scale coronal field is dominated by 
the polar fields. 



4. Coronal Model with Anisotropic Gas Pressure 

At cycle minimum, the corona presents a relatively ordered structure with high latitude 
coronal holes and an equatorial streamer belt. The corona is essentially axisymmetric and 
stable for many months (e.g., Gibson 2001). In the following we describe an axisymmetric 
model of the corona; an earlier version of this model was described by Vasquez, Raymond 
& van Ballegooijen (1999) and Vasquez, van Ballegooijen & Raymond (1999). We use a 
spherical coordinate system (r, 6, 0) and assume that all scalar quantities are independent 
of azimuth <p. The magnetic field B is assumed to lie in the meridional plane: 

„ > I fldA^ dA\ 
B = V x A = — - -— e r - — e e , (6) 
rsmtf \r ot> or J 

where A = A^e^ is the vector potential, and A = A^rsinO is the field-line variable (note 
that A = constant along field lines). This variable increases from A = along the polar axis 
to A = A eq at the equator (r = R Q , 9 = ir/2). There is a critical field line A c that forms the 
boundary between open and closed magnetic regions; the field lines with < A < A c are 
open, while those with A c < A < A eq are closed. Note that A c /A eq is the fraction of surface 
flux that is open. In this paper we assume A c /A eq = 0.80. The magnetic configuration is 
illustrated in Figure 4. 

Rather than solving an energy equation for the coronal plasma, our semi-empirical model 
is based on observed temperature distributions. The electron and proton temperatures (T e , 
T p \\ and T p j_) are interpolated between the pole and equator: 

T(r, A) = T pole (r) + [T equa (r) - T poie (r)]$(A), (7) 

where T po ; e (r) and T equa (r) are the observed profiles described in §3, and $(A) is a function 
of the field-line variable that increases monotonically from $(0) ^ at the pole to &(A C ) 1 
at the equator. We use the following expression for <&(A): 

HA) = l + e W [-(A-A h )/A w y (8) 
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where Ah is the mid-point of the transition = |], and A w measures the width of the 

transition. In this paper we use Ah/A c = 0.70 and A W /(A C — Ah) = 0.15. The function &(A) 
is shown in Figure 5. 

The coronal magnetic and velocity fields are computed by solving the momentum equa- 
tion: | 

V • (P + pvv) = -pV$ g + — (V x B) x B, (9) 

where P is the pressure tensor, p is the mass density, v is the outflow velocity, 3> 9 (r) = 
—GM & /r is the gravitational potential. The steady-flow condition requires that v || B, and 
conservation of mass requires that pv/B is constant along field lines. The pressure tensor is 
anisotropic: 

P=P||SS+p ± (l-ss), (10) 

where I is the unit tensor, s is the unit vector along B, and p\\ and p± are the parallel and 
perpendicular pressures: 

PH = N e k B {T e + T p]] ), (11) 
P± = N e k B (T e + T p± ). (12) 

Here we assume a pure hydrogen plasma. The left-hand side of equation (9) can be written 

as 

V.(P + pvv) = Vp, + ^(!)s + *4 (13) 

where s measures distance along the field lines, B = |B|, 5p = p\\ — p± + pv 2 , and we used 
V • B = 0. The component of the momentum equation (9) parallel to s reads 

dp± +B d(6p\ = _ p ^ (14) 



ds ds \B J ds ' 

and the solution of this equation will be further discussed in §4.2. We now introduce the 
function P±(r, A) describing the dependence of p± on the field-line variable: p±(r, 6) = 
P±[r, A(r, 9)}. Then the gradient of p± is 

where the partial derivatives on the RHS are taken at constant A and constant r, respectively. 
It follows that dp±/ds = (dP±/dr) cos7, where 7 is the angle between s and e r , and similarly 
d&g/ds = (d&g/dr) cos 7. Inserting these expressions into equation (14), we obtain 

dP ± d§ q B d (5p\ , , 

or dr cos 7 as \ B J 
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and combining equation (9), (13), (15) and (16) yields 

&Pi , B d[Sp\^ _ . ds (VxBL _ _ , 

-VA + — -| (scos 7 - e r ) + Sp— = ^— J -^(-B e e r + B r e e ). (17) 



dA cos 'j ds \ B J ds Att 

In streamers the pressure anisotropy and outflow velocity are small, so the terms involving 
Sp can be neglected. The same is true in coronal holes because in these regions the magnetic 
field is nearly radial and the vectors (s cos 7— e r ) and ds / ds vanish for radial fields. Therefore, 
Sp plays only a minor role in the perpendicular force balance, and we will neglect its effects 
in equation (17). Then the perpendicular force balance (17) reduces to 

(V x B) . dP ± 

-^^ = ™ne—, (18) 

where we used equation (6). The solution of equation (18) will be discussed in §4.1. 

We use an iterative method for solving equations (16) and (18). The method can be 
summarized as follows. We first construct a preliminary pressure model P±(r,A) by inter- 
polating the observed perpendicular pressure between pole and equator, using an expression 
similar to equation (7). Then we compute the field-line variable A(r, 9) by solving the per- 
pendicular force balance as described below in §4.1; this gives us a first guess for the shape 
of the magnetic field lines. Then we solve the parallel force balance separately along many 
field lines as described in §4.2. This yields the mass density p(s) and velocity v(s) along 
each field line, which are then remapped to produce the density p(r,9) and velocity v(r,9). 
We also update the perpendicular pressure P±(r,A), which is then used in the next itera- 
tion to recompute the field-line variable A(r,9). This process is repeated until convergence 
is achieved. Note that in the first iteration P±(r,A) is specified analytically, whereas in 
subsequent iterations P±(r, A) is obtained by numerical interpolation. 



4.1. Perpendicular force balance 



To solve equation (18) for the perpendicular force balance, we introduce the following 
integral over the coronal volume: 



/•t/2 




/o 


JRq 



52 

8^ 



-P±(r,A) 



r 2 sin 9 dr d9, 



(19) 



where only one hemisphere of the Sun is considered, and R m ax = 10 Rq is the outer radius 
of the computational domain. The magnetic field B is related to A^ via equation (6), and 
A = A^r sin 9. At each step in our iterative procedure, P±(r, A) is a known function, so 
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L depends only on the unknown function A^(r, 9), or equivalently, A(r,9). We assume the 
following boundary conditions for A(r, 6): 

A(R Q ,6) = ^(l-cos^fl), (20) 
p + 1 

OA 

^(Rmax,0) = 0, (21) 

or 

A(r,0) = 0, (22) 

A(r,n/2) = A c for r > r cusp , (23) 
dA 

— (r,7r/2) = for r < r cusp . (24) 

Equation (20) is derived from the surface flux distribution, equation (5), and the magnetic 
field is assumed to be radial at the outer boundary, r = R max . Equations (23) and (24) 
describe the boundary conditions at the equator. For r > r cusp , there is a current sheet at 
the equator separating the open fields from the two hemispheres, hence the magnetic field 
just above the equator is radial. For r < r cusp the magnetic field is closed over the equator, 
so the field is perpendicular to the equatorial plane. 

We now show that the function A^(r, 6) for which L reaches its minimum value is 
a solution of equation (18). Let (L4^(r, #) be an arbitrary variation of A^,{r,0), then the 
corresponding change in L is 



f7r/2 r R„ 

5L 



/■7T/2 ptimax 
JO J Rp, 



(VxB)^ . dP± 
— rsmd- 



4tt OA 



5 A^r, 9) r 2 sin9 dr d9, (25) 



where the boundary conditions (20) - (24) have been used. At the minimum, 5L vanishes 
for any function 8A^{r,0), so the quantity in square brackets in equation (25) must vanish. 
This is precisely the condition for perpendicular force balance, equation (18). 

The function A^r, 6) is discretized on a non-uniform grid in r and 6 with 180 points in 
each direction. The minimization of L is performed using the Polak-Ribiere version of the 
conjugate gradient method (Press, et al. 1992). This is an iterative method for adjusting 
the values of A^ at the 180 2 grid points until the minimum of L is reached. The cusp radius 
r cusp is also allowed to change in order to obtain the lowest possible L, but the amount of 
open magnetic flux (A c /A eq ) is held fixed. 
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4.2. Parallel force balance 

The parallel component of the momentum equation (14) can be written in the following 



form: 



ds 



p , 9 ,,i 9 1 dB d& a , , 

|H + - 2 )J+pci-— = -p^, (26) 



B ds ds ' 

where cm and cj_ are the parallel and perpendicular sound speeds: 



„2 _ P||,-L _ k B(Tp \\,±+T e ) 



c 2 nu± = ^^p^ • (27) 
p m p 

Using pv/B = constant along field lines, we can eliminate the density from equation (26): 
where 

f(s)~ GM& <?(*) [ Scl{s ' )dB d S > (29) 

f{s) = 75T " c " (s) " y rfs • (29) 

At each step of our iterative procedure we have estimates for B(s), r(s), T p \\(s), T p ±(s) and 
T e (s) on any open field line, which allows us to compute the function f(s). According to 
equation (28), the sonic point (v = cy) is located at an extremum of f(s). The velocity v(s) 
is found by inward and outward integration of equation (28), starting at the sonic point. 
In practice we find that, in order to obtain a valid solution over the entire height range 
(Rq < r < Rjnax), the sonic point must be located at the global minimum of f(s). The 
density p(s) is computed from mass flux conservation and the boundary condition for the 
density at the coronal base (see §3). This process is repeated for 180 different open field 
lines, and the results are remapped to obtain the density p(r,9) and velocity v(r,0) on the 
(r, 9) grid. 



5. Results 

We iterated the parallel and perpendicular force balance equations as described in §4. 
The process was repeated until the maximum change in P±(r, A) between iterations is less 
than 1%; this required 20 iterations. The results for the first and last iterations are shown 
in Figures 6a and 6b, respectively. The thin curves represent magnetic field lines, and the 
thick curves are contours of (3(r, 9), the ratio of perpendicular gas pressure p±_ and magnetic 
pressure B 2 /8n. The cusp height in the final solution is r cusp 1.95i? . Note that the 
magnetic structure changes little between the first and last iterations. 
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Our model yields (3 <C 1 in the polar regions, but (3 > 1 in the equatorial streamer, 
especially near the streamer cusp. This supports the idea that the streamer is magnetically 
contained by the strong polar fields that surround it on either side (Suess, Gary & Nerney 
1999). We find that f3 > 1 throughout the closed-field region of the streamer. However, the 
variation with height is not monotonic: (3 has peaks at both the cusp and the streamer base, 
and lower values at intermediate heights. Our results are similar to those of Li et al. (1998), 
who estimated (3 based on UVCS/ SOHO and SXT/Yohkoh observations of streamers (July 
1996) in combination with potential field extrapolation of the photospheric magnetic field. 
Their estimates indicate (3 ~ 5 at 1.15 R Q and (3 ~ 3 at 1.50 R Q) similar to the values found 
here. High values of (3 were also found in MHD models that include heat and momentum 
deposition in the corona (Wang A. H. et al. 1998a; Suess et al. 1996). 

Figure 7 shows the outflow velocity for the final model. The bright region is the fast solar 
wind emanating from the coronal hole, and the dark region is the slow wind that flows along 
the open field lines within the streamer. The maximum velocities are about 450 km/s for 
the fast wind and 190 km/s for the slow wind. This is somewhat smaller than the observed 
in-situ values at 1 AU (800 and 400 km/s, respectively). We attribute this difference to the 
fact that our model does not include any momentum deposition effects other than the mirror 
force. However, the ratio of fast and slow wind speeds is about 2.5 in our model, consistent 
with in-situ measurements. 

Figure 8 shows three quantities measured across magnetic field lines: the temperature 
anisotropy at 3 R & , the radial position r s of the sonic point, and the asymptotic wind speed. 
Note that the variation of wind speed closely follows that of the temperature anisotropy 
(compare top and bottom panels). This indicates that the fast wind is driven by the high 
perpendicular temperatures in the coronal hole, and the decrease in wind speed at the edge 
of the hole is mainly due to the decrease in temperature anisotropy. The middle panel 
shows that there is a sudden jump in sonic-point radius once the transition from high to low 
temperature anisotropy is nearly complete. The jump occurs at A/A c = 0.77 and is due to 
the appearance of a second (lower) minimum in the function f(s) at a height well beyond 
the streamer cusp (see below). As a result, the sonic-point radius changes discontinuously 
from about 1.5 R Q in the coronal hole to about 5.5 R & in the streamer. Such large values for 
the sonic-point height in the streamer were found earlier by Wang Y.-M. (1994) and Chen 
& Hu (2001). Most of the decrease in wind speed occurs for A/A c between 0.6 and 0.77, 
well before the jump in sonic-point height, and the asymptotic wind speed does not change 
significantly at the jump. 

To show the transition between fast and slow wind more clearly, Figure 9 shows various 
quantities along field lines near the fast-slow boundary (the same line styles are used in all 
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panels). Panel (a) show the shapes of the field lines in the meridional plane. Panel (b) 
shows the function f(s) defined in equation (29) and plotted as function of radial distance 
r for four field lines. Note that the triple- dot- dashed curve has a minimum at about 1.5 R Q , 
whereas the full curve has two minima, one at 1.7 R & and another at 5.5 R & . As shown in 
§4.2, these minima indicate possible positions of the sonic point (v = c\\). Panel (d) shows 
the radial positions r m i n of such minima as function of field-line variable A. In general, a 
global solution of the wind equation (28) can be found only when the integration is started 
at the global minimum, i.e., the minimum with the lowest value of /(r min ). Therefore, as we 
move from the triple- dot- dashed field lines to the dot-dashed line, there is a discontinuous 
jump in the height of the sonic point. 

Figure 9c shows the outflow velocity v(r) along field lines near the fast-slow boundary. 
The dotted curve in this panel is the wind solution along the last open field line, A = A c . 
As we move from the coronal hole into the streamer, the wind speed at large height closely 
follows the decrease in perpendicular temperature, which mainly occurs between the two 
full-thick profiles. However, a further decrease occurs at lower heights (r ~ 2 R Q ) between 
the triple- dot- dashed and dot-dashed field lines due to the jump in height of the sonic point. 
We conclude that, unlike the terminal speed, the low-height behavior of v(r) is strongly 
affected by the height of the sonic point. Such a region of stagnated flow was found earlier 
by many authors (e.g., Wang Y.-M. 1994; Suess & Nerney 1999, 2002; Chen & Hu 2001). 

Figure 9e shows the non-radial expansion factor of the slow- wind flow tubes, f ex p{ r ) — 
(R & /r) 2 B(R Q )/B(r), where B(R Q ) in the field strength at the coronal base for each field 
line. Note that / exp (r) is a non-monotonic function of radius for these field lines. The peak 
in / exp (r) is due to the fact that these field lines pass close to the streamer cusp where B xs 0. 
In contrast, f cxp ( r ) increases monotonically for field lines in the coronal hole (not shown). 
It is remarkable to see that this "cusp effect" occurs at large distance from the cusp: the 
peak in / exp (r) occurs for all field lines within about 1 R & from the equatorial plane (see 
Figure 9a). The triple- dot- dashed field line at which the peak in / exp (r) first develops is close 
to the field line where the jump in sonic-point height occurs. This is consistent with the 
suggestion by Suess & Nerney (1999, 2002), Chen & Hu (2002) and others that the existence 
of the slow wind is due to the non-monotonic behavior of / exp ( r )- We suggest here that the 
decrease in perpendicular temperature as we approach the equatorial plane also plays an 
important role. 

To determine whether the "cusp effect" is due to the effects of gas pressure on the 
expansion of the flow tubes, we computed a partially open potential magnetic field with the 
same photospheric boundary conditions [see equation (5)] and a thin current sheet in the 
equatorial plane for r > 1.95_R . The calculation is based on the work of Low (1986). First, 
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the surface flux distribution is expressed in terms of four Legendre polynomials P n (x) with 
n — 1, 3, 5 and 7. Then the field is extrapolated into the corona using equations (B3) and 
(B4) of Low (1986) for n = 1 and 3 (for n = 5 and 7 the effect of the current sheet can 
be neglected). We found that f ex p{ r ) has a peak for all open field lines that pass below the 
point (x, z) = (3.0,1.4) R & , where x is the distance from the rotation axis, and z is the height 
above the equatorial plane. Therefore, the effect of the streamer cusp -located at the point 
(x,z) = (1.95,0) R Q - is present even in a potential field, and its spatial extent is similar to 
what was found in the numerical model (i.e., the effect is not due to the finite plasma (5). 

For the value of Ah chosen in this paper (Ah/A c = 0.7), the transition between fast 
and slow wind far from the Sun occurs at latitudes of ±20°. This is consistent with Ulysses 
observations taken at the time of cycle minimum (Suess et al. 1999). At lower heights, the 
low speeds in the streamer legs found in our model are consistent with UVCS measurements 
of 5+ outflow velocity using the Doppler dimming technique (Strachan et al. 2002). In 
a latitudinal scan at 2.33 Rq, Strachan et al. (2002) found no measurable outflow velocity 
in the streamer (v < 20 km/s), but a steep rise in velocity occurs just beyond the bright 
streamer legs. The observed latitudinal half-width of the streamer at 2.33 R & , as defined by 
the points where the outflow velocity equals 100 km/s, is about 20°. Our model predicts a 
half-width of 22°. Therefore, the predicted width of the streamer is roughly consistent with 
observations at both large and small heights. 

It should be noted that we use a realistic photospheric flux distribution that is peaked 
at the pole, B r (R & ,9) = 10 cos 7 9 [G], whereas Suess, Gary & Nerney (1999) and Chen & 
Hu (2001) assume a bipolar distribution, B r (R & ,9) oc cos 9. For a given fraction of open 
magnetic flux, the assumed flux distribution has a significant effect on the size of the polar 
coronal hole, the latitudinal width of the streamer core, and the magnitude of the non-radial 
expansion factor. However, potential field models indicate that the spatial extent of the cusp 
effect is similar in the two cases. 



6. Comparison with Visible Light and Lya Observations 

The coronal density model can be used to predict the polarization brightness (pB) of 
the visible light that is scattered by free electrons in the corona (Thomson scattering). The 
pB is an integral along the LOS (van de Hulst 1950): 

pBM s bJ NAx) (^f (l^M^M dx , ,30) 

16 J L os V r ) 1 - u'/3 

where x is the distance along the LOS, r is the projected radial distance from Sun center, 
r = ^/r 2 + x 2 is the true radial distance, a T = 6.65 x 10~ 25 cm 2 is the Thomson scattering 
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cross section, B Q = 1.97 x 10 10 erg cm~ 2 s _1 sr _1 is the mean disk intensity, and u' = 0.63 
is the limb darkening coefficient. The functions A(r) and B(r) are given by van de Hulst 
(1950). We used our axisymmetric model of the electron density distribution to compute 
synthetic pB images for various angles a between the LOS and the solar equatorial plane. 
Figure 10 shows the results for a = and a = 30°; the latter represents a view of the Sun 
from above the ecliptic plane. We have over-plotted selected field lines in the plane of the 
sky (white curves) and contours of the ratio of pB and its maximum value at the equatorial 
base (black curves). The equatorial streamer clearly stands out in the image with a — 0, but 
is less distinct for a = 30°. 

To compare these results with observations, Figure 11a shows the radial pB profiles 
along the pole and equator, and Figure lib shows latitudinal profiles at three different 
heights (r = 1.15, 1.5 and 2.5 R Q ). The full curves correspond to the edge-on view (a = 0). 
The symbols represent pB measurements obtained with the Mauna Loa Mark 3 coronagraph 
and LASCO / SOHO (Gibson et al. 1999; Guhathakurta et al. 1999) and with the visible-light 
channel on UVCS/ SOHO (Cranmer et al. 1999). The data were obtained in 1996 July and 
August when the equatorial streamer belt is seen approximately edge-on (Raymond et al. 
1997). Our model fits the observations quite well, although the predicted latitudinal width 
of the streamer at 2.5 R & is somewhat larger than observed. The dashed curve in Figure 11a 
shows the radial profile along the pole for the out-of-ecliptic view (a = 30°). Note that there 
is a large increase in pB compared to the edge-on view. This is due to the fact that the LOS 
crosses the enhanced density equatorial streamer, and is not representative of the available 
observations of the polar pB profile. 

We also computed the intensity of the strongest UV coronal emission line, H I 1216 A 
(Lycc), which is formed almost entirely by resonant scattering of chromospheric Lja radiation 
by neutral hydrogen atoms in the corona (Noci, Kohl & Withbroe 1987; Cranmer et al. 1999). 
The integrated intensity of the scattered Lya is 

I = ^B 12 [ N H (x)Q(r) [ I & (v + 5v)<j>{y - u ) dv dx, (31) 
47r Aos J 

where x measures distance along the LOS, Nh(x) is the neutral hydrogen density, r is the 
radial distance from Sun center at any point along the LOS, tt(r) = [1 — a/1 — (R Q /r) 2 ]/2 is 
a geometric dilution factor, i©(z^) is the chromospheric intensity (in erg cm _2 s _1 Hz _1 sr _1 ) as 
function of frequency v (in Hz), uq = 2.47 x 10 15 Hz is the line center frequency, bv = vqv t /c 
is the Doppler shift, v r is the radial component of the outflow velocity, 4>(v — u ) is the 
scattering profile of the coronal atoms (due to their thermal and non-thermal velocities), 
and B\2 is the Einstein coefficient. For a derivation of this expression and discussion of the 
relevant approximations, see Noci, Kohl & Withbroe (1987). The neutral hydrogen density 
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N H is computed using the collisional ionization rates of Scholz & Walters (1991) and the 
recombination rates of Hummer (1994). 

We computed Lja images for two different viewing angles, using densities, temperatures 
and velocities from the coronal model. The results are shown in Figure 12. Although these 
images look similar to those for polarization brightness, the Lja intensity is somewhat 
sensitive to the outflow velocity. This is shown more clearly in Figure 13 where we plot the 
radial profiles of Lja intensity with and without the Doppler dimming effect (full and dashed 
curves, respectively). Note that there is a significant difference between predicted intensities 
with and without dimming for the polar region, but not for the equatorial region. Therefore, 
unlike pB measurements Lya intensities provide constraints on the outflow velocity in the 
acceleration region of the solar wind. 

Figure 13 also compares our model predictions with measurements from UVCS / SOHO. 
The diamonds represent measurements taken along the streamer axis, based on our reanalysis 
of data obtained on 1996 July 26 (Raymond et al. 1997). The triangles are measurements 
along the polar axis based on observations from November 1996 to April 1997 (Cranmer et 
al. 1999; Dobrzycka et al. 1999). Note that the predicted Lya intensities fit the observations 
quite well. This suggests that the model gives a reasonable representation of the velocity 
field in the inner part of the corona where the main acceleration of the solar wind takes 
place. 

7. Discussion 

We developed a stationary, axisymmetric MHD model for the global corona. The model 
includes a description of the temperature, density, velocity and magnetic field as function of 
latitude and radius up to 10 Rq from sun center. The velocity and magnetic fields are ob- 
tained by solving the parallel and perpendicular force balance equations, including the effects 
of inertia, anisotropic gas pressure, gravity and Lorentz forces. The temperature models are 
based on observational data from UVCS/ SOHO, and the magnetic flux distribution at the 
coronal base is taken from NSO/Kitt Peak synoptic maps. The model reproduces the main 
features of the global corona at the time of cycle minimum. We find that the high perpen- 
dicular temperature of the protons in the coronal hole plays a major role in driving the fast 
solar wind. In the streamer we find low outflow velocity and high plasma (3, consistent with 
earlier results (Suess et al. 1996; Wang A. H. et al. 1998a; Li et al. 1998). 

In our model the wind equation (28) is solved separately for many field lines. The sonic 
point along each field line occurs at a minimum of the function f(s) that appears on the 
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RHS of the wind equation [see equation (29)]. The transition from fast to slow wind occurs 
at an open field line characterized by A h = 0.7 A c , where A c is the boundary between open 
and closed fields. The value of was adjusted to obtain the correct latitudinal width of 
the slow-wind region both at large distance from the Sun (Suess et al. 1999) and in the 
corona at 2.33 R® (Strachan et al. 2002). In our model this transition is associated with 
two effects: (1) the decrease of proton perpendicular temperature T p ±(r, A) as we approach 
the equatorial plane, and (2) the appearance of a peak in the non-radial expansion factor 
fexp(f) fc> r held lines that pass close to the streamer cusp ("cusp effect"). The cusp effect 
is present at surprisingly large distances from the cusp (~ 1 R Q ) and is present even in 
potential-field models, so it is not a consequence of finite plasma (5. The combination of 
decreasing temperature anisotropy and cusp effect causes the global minimum of f(s) (and 
therefore the sonic point) to occur well beyond the cusp, and produces low outflow velocity 
near the cusp. These results are consistent with the suggestion by Noci et al. (1997) that the 
slow wind is due to special properties of the geometric spreading along the open field lines 
that pass near the streamer core (also see Wang & Sheeley 1990; Wang Y.-M. 1994; Wang 
Y.-M. et al. 1997b; Chen & Hu 2001, 2002; Suess & Nerney 1999, 2002). 

Our model does not provide a physical explanation for the temperature decrease at the 
fast-slow boundary, and therefore cannot explain why the boundary occurs at Ah = 0.7 A c . 
Understanding the physics of the fast-slow transition will require more detailed analysis of 
the energy balance of the coronal plasma, including the physical processes by which the tem- 
perature anisotropy of the protons is maintained. We speculate that proton perpendicular 
heating (by dissipation of transverse MHD waves) occurs in both the fast and slow winds, 
perhaps at roughly equal rates. However, the resulting temperature anisotropy T p ±/T p \\ may 
be quite different in the two cases. In the fast wind, proton-proton collisions are less fre- 
quent due to the lower density, so the deviations from Maxwellian velocity distributions are 
larger than in the slow wind. Clearly, to understand why the fast-slow transition occurs at 
Ah = 0.7A C will require multi-dimensional models of wave heating and energy balance such 
a those developed by Chen & Hu (2001). 
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Fig. 1. — Electron density profiles derived from Skylab observations for the pole (thin curve) 
and the equator (thick curve). 

Fig. 2. — Temperature profiles derived from coronal observations, for the pole and the 
equator. 

Fig. 3. — Photospheric magnetic field as function of latitude at the time of the last cycle 
minimum (July 1996). The diamonds are measurements from Kitt Peak synoptic maps. The 
full curve is a fit to the data. 

Fig. 4. — Sketch of the dipolar configuration. 

Fig. 5. — The function used for interpolating temperatures between pole and equator. 

Fig. 6. — Magnetic field and plasma pressure for the first iteration (top) and for the final 
model (bottom). The thin curves are magnetic field lines (contours of A), and the thick curves 
are contours of f3 = 8np±/B 2 . 

Fig. 7. — Outflow velocity in the meridional plane. The brightness level increases with 
magnitude of the outflow speed, v, as shown by the scale on the left side of the figure (v 
in km/s). The solar rotation axis is along the left axis of the plot, and the solar surface is 
indicated by the thick semi-circle. The white curves are selected field lines labeled with their 
value of A/A c . 

Fig. 8. — Various wind properties as function of field-line variable A. Top panel: temperature 
anisotropy (T p ±+T e ) / (T p \\+T e ) at r = 3R Q . Middle panel: radius of sonic point where v = c\\ 
(in R Q ). Bottom panel: wind speed at R max = 10 Rq (in km/s). The thin solid vertical line 
in each panel indicates the field line A^ where the latitudinal temperature gradient is largest, 
and the thick vertical lines indicate the width of the transition (A h ± A w ). The dot-dashed 
vertical line indicates the field line at which there is a sudden transition in the radius of the 
sonic point. 

Fig. 9. — Transition between fast and slow wind, (a) Shape of various field lines in the 
meridional plane, (b) Functions f(r) that appear on the RHS of the wind equation (28) for 
different field lines, (c) Outflow velocity v(r) along field lines. The dotted curve corresponds 
to the last open field line, A = A c . (d) Radial positions of the mimima of f(r) as function 
of field-line variable A (diamonds). The vertical lines show the A- values of the field lines 
shown in other panels, (e) Non-radial expansion factors f C x P (r) for different field lines. The 
same line styles are used in all panels. 
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Fig. 10. — Images of visible light polarization brightness (pB) as predicted from the coronal 
model. Black curves are contours of \og(pB/pB max ), where pB max is the value of pB at the 
coronal base. The white curves are selected field lines in the plane of the sky (labeled by 
A/ A^. The dark semi-circle represents the solar disk. 

Fig. 11. — Left panel: Visible-light polarization brightness of the corona, pB/B Q , as function 
of radial distance from sun center. The solid curves show the predicted equatorial (thick) 
and polar (thin) profiles for an edge-on view. The dashed curve is the predicted polar profile 
for a view from above the ecliptic plane (a = 30°). The symbols are measurements for 
the pole (triangles) and equator (diamonds) from Gibson et al. (1999), Guhathakurta et al. 
(1999), and Cranmer et al. (1999). Right panel: Predicted polarization brightness as function 
of co-latitude at projected radii of 1.15, 1.5, and 2.5 R®. The symbols are corresponding 
measurements from the above-cited references. 

Fig. 12. — Images of Lya intensity / as predicted from the coronal model. The black curves 
are contours of log(J / I max ) , where I max is the maximum value of / at the coronal base. The 
white curves are selected field lines in the plane of the sky (labeled with A/A c ). The dark 
semi-circle represents the solar disk. 

Fig. 13. — Radial profiles of Lja intensity I. The solid curves are the predicted intensity 
along the pole (thin solid curve) and along the equator (thick solid curve) for the edge-on 
view (a = 0). The dashed curves are the predicted profiles without Doppler dimming [i.e., 
with v r = in equation (31)]. The triangles and diamonds are intensity measurements from 
UVCS/ SOHO for the pole and equator, respectively (see text). 
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Table 1. Numerical values of the parameters for the temperature models. 





T p ± 


T p \\,T e 






at Pole 


at Pole 


at Equator 




8 x 10 5 K 


8 x 10 5 K 


8 x 10 5 K 


Ti 


3 x 10 6 K 


0. 


0. 


a 


0. 


0. 


0.1 


b 


0.23 


0.47 


0.33 


a 


0.7 


0.7 


0.55 




6.6 


6.6 


6.6 


Ar 


1R Q 
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